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TXTLE OF TEE X»VEKTI0N: 

TEDS KSSIHEERXHG OF BX0SCTROST&TXC IBTEE&CrXOHS A$ 
MET&L XOJ* BINDING SITES FOR THE 
STABX LT.2ATXOJ* OF PROTEINS 

This is a continuation-in-part application of tUS, 
serial Ho. 034,965, filed April 6, 1937, th* contents of 
which are fully incorporated herein by reference. 

\hm present inyantion relates to designing proteins 
with enhanced stability through the alteration of the 
proteins* divalent metal ion binding sites. 

Background o f the x .nv-ent.i.o.n 

Proteins are linear polymers of amino acids. Since 
the polymer i nation reaction which produces a protein 
results in the lose of one molecule of water from each 
amino acid, proteins are often said to bo composed of 
amino acid "residues." Natural protein molecules roay 
contain as many as 30 different types of amino acid 
residues, each of which contains a distinctive side 
chain. The sequence of axaino acids in a protein defines 
the primary structure of the protein, 

Proteins fold into a three-dimensional structure. 
The folding is determined by the sequence of amino acids 
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and by the protein's environment. remarkable 
properties of proteins depend directly from the 
protein ! s three-dimensional conformation. Thus, this 
conformation determines the activity or stability of 
enay^es, the capacity and specificity of binding 
proteins, ana the structural attributes of receptor 
molecules. Because the three-dimensional structure of a 
protein molecule is- so significant., it has long been 
recognised . that a means for stabilising a protein's 
three-dimensional structure would he highly desirable. 

The three-dimensional structure of a protein may tee 
determined in a number of ways. Perhaps the best known 
way of determining protein structure involves the use ©f 
the technique of • 3t~ray crystal iography* An excellent 
general review of this technique can be found in 
^Y^jc^l.^toMi^ry, van Hold*, K..E-. (Prentice-Hall, 
KJ (1971) pp221~23S} which reference is herein 
incorporated by reference. Using this technique, it is 
possible to elucidate three-dimensional, structure with 
remarkable precision. It is also possible to probe the 
three-dimensional structure of a protein using circular 
diuhroism, "light scattering, or by measuring the 
absorption and emission of radiant energy (Van Holde, 

* - ^" , NJ (1971) } , 
Additionally,, protein structure may he determined 
through the use of the t< neutron de fraction,, 

cr by nuclear magnet- ic resonance { Ehysicai C hem istry, 

4th !d« Moore, W.J.., Prentice-Hall, HJ (1972) which 
reference is hereby incorporated by reference) . 

The examination Of the 'tbree-dimens i one I structure 
of numerous natural proteins has revealed a number of 
recurring patterns. Alpha helices, pa rails! beta 
sheets, and anti-parallel beta sheets are the most 
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common patterns observed, An excellent description of 
such protein patterns is provided by Dicker son, R.S. , at 
al* In; -The structure and Action of Frpte,fiis f 
Benjamin, Inc., CA (1969) * The assignment; of each amino 
acid to one of these patterns defines the secondary 
structure of the protein. The helices, sheets and 
turns of a protein's secondary structure pack together 
to produce the three-dimensional structure of the 
protein. The three-dimensional structure of ssany 
proteins cay be characterised as having internal 
surfaces (directed away from the aqueous environment in 
which the protein is normally found) and external 
surfaces (which are in close proximity to the aqueous 
environment) . Through the study of many natural 
proteins, researchers have discovered that hydrophobic 
residues (such as tryptophan, phenylalanine, tyrosine > 
leucine, isoieueine, valine, or methionine) ere most 
frequently found on the internal surface of protein 
molecules. In contrast,, hydrophilic residues (such as 
asperate, asparagine, glutamata, giutamins, lysine, 
erginine f histidine, serins, threonine, glycine, and 
proline) are most frequently found on the external 
protein surface. The amino acids alanine, glycine,, 
serine and threonine are encountered with equal 
frequency on both the internal and external protein 
surfaces. 

Proteins exist in a dynamic equilibrium between a 
folded, ordered state and an unfolded, disordered state. 
This equilibrium in part reflects the short range inter- 
actions between the different segments of the poly- 
peptide chain which tend to stabilise the protein's 
structure, and, on the other hand, those thermodynamic 
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forces which tend to promote the random! g at ion of the 
molecule. 

Metal ions have long been known to stabilise 
proteins -fey binding at specific sites in the tertiary 
structure. Indeed* ©any proteins isolated from therm©- * 
philie organisms been shown to contain calcium ion 

binding sites. These sites are responsible for the 
enhanced stability that these proteins require in order 
to function in the elevated temperature milieu in which 
they are found. For example, thermolysin, a neutral 
protease from the thermophilic organism Bacillus 
thermgprc jiyj - /as found to contain fowr calcium 
binding sites. By studying the calcium ion dependence 
for the rate of thermal inachivation*. it has been 
possible to estimate the collective contribution of 
these sites to the G of unfolding to be between 8.1 and 

9.2 Kcal/mole (vbrrdouw et el, , Biochemistry 15:3716- 

3724 (1976)}. 

For example, Serparsu e-t„ai. f | -cry 26: 12S3~ 

1300 {1987} replaced the carhoxylata ligand of residue 
aspartic acid 40 in the calcium binding site of 
staphylococcal nuclease with a glycine resulting in a 
7.4 fold decrease in affinity for calcium ions. 
Similarly , Tsuju g&j&L* > Sms^M^LJm^t (HSU W&> 
83 107-8 XXX {1986} changed glycines in the vicinity of 
the calcium binding loops in the photoprotein aeguorin 
to positively charged arginines to cause a decrease in 
th® affinity of metal binding at these sites. It would, 
however, be desirable to design proteins with increased 
stability by increasing the binding affinity of divalent 
metal cations at a specific site in the protein through * 
specific' amino acid substitutions. 
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The present invention provides a method for 
redesigning proteins to incs**** the stability of the 
protein fey altering amino acid residue fs) that ere in 
close proximity to the protein's metal ion binding 
site(s). In detail, the method of this invention 
involves the following steps : 

1, correlating the metal ion binding structural 
site with some parameter of protein stability? and 

2, Altering the metal ion binding site through 
the substitution, insertion,, or delation of amino acid 
residues in close proximity to that site so that the 
electrostatic attraction interaction between the amino 
acids end the metal ion is increased. 

in identifying the metal ion binding site for 
correlation with a parameter of protein stability, one 
may analyse the three dimensional structure cf a protein 
of interest or some evolutionalXy related variant* The 
three dimensional structure may he obtained from 
published sources, or determined by Known methods of X- 
ray crystal Xography * 

In altering the -metal ion binding site through the 
substitution, insertion., or deletion, of amino acid 
residues, the following steps are considered in the 
method of this invention: 

(a) Selecting amino acids for substitution., 
insertion, or deletion as close as possible to the metal 
ion binding site, but without introducing steric 
hindrance ? 

(b) Preferentially selecting for substitution., 
insertion, or deletion those assino acids that are not 
conserved in evolutional Xy related homologous proteins? 
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(c) changing positively efcarged amino acid 
residues that are close to the bound metal ion. to 
uncharged or negatively charged aaino acid residues? and 

(d) Preferentially using computer-assisted 
graphics to simulate the introduction of the new amine 
acid residues. 

Farther, the invention is directed to proteins that 
have been redesigned to have enhanced stability 
according to the methods- of this invention. 

SrJLejg t>eggritttia a-ffif„tfefe.. Figures 

E&gMSMJk shows a i» 3i resolution X-ray structure of 
subtil i sin BPU 1 in stereo. The two calcium binding 
sites are labeled CaA and CaB> 

Fiqrare ,3 shows -detailed structural infdotati.cn for 
the calcium A {CaA} binding site in stereo, 

ligaiaJ, shows detailed structural information for 
the calcium 8 (CaB) binding site in sterag, 

EiSSS&dl shows calcium ion d«pend«nc« for the rate 
of thermal inactivation of subt il lain at 65 5 c. The data 
for subtilisin 7172 is denoted by {*) and for the wild- 
type subtilisin BPB* by («>.. The < * } denotes the 
effect of Hgci 2 on the subtilisin BPS?* » The solid linos 
represent theoretical titration curves for single metal 
ion binding sites having the affinity (%} shown, 

shows calcium ion dependence for the rata 
of thermal inactivation of subtiiisin 7148 at 65 'C 

Definitions 

The following definitions are used in describing 
the invention, 

ffrgta in 
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A heteropolyaer mad® by living cells and composed 
of anino acids, h typical protein comprises 100 to 1000 
amino acids. The ex-act sequence of amino acids deter- 
mines the structure and faction of the protein. 

one of the twenty naturally occurring compounds 
which are the building blocks of proteins* The natural 
amino acids are usually abbreviated to either three 
letters or one letter according to Table l. The amino 
acids are joined head to tail to form a long main chain, 
Each kind of amino acid has a different side group. 
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Table 1. Amino acid names and abbreviations. 



Amino Acid Three latter cod® Single letter code 



Alanine Ala 







1 Pi tlXO id 


Asp 




Asn 










~r ine 


Gin 


Glycine 


Sly 


Histidine 


Mis 


Xso leucine 


lie 


Leucine 


Leu 


Methionine 


Met 


Phenyls 1 an ine 


Phe 


Proline 


Pro 


serine 


s«sr 


Threonine 


Thr 


Try p t oph a n e 


Trp 


Tyrosine 


T$x 


Valine 


Val 



Ail amino acids &ave the same atoms in the stain 
chain and differ only in the side chains, The main- 
chain atoms are a nitrogen, two cartoons, and one oxygen. 
The . first atom is the nitrogen, called simply if. The 
next ..atom is a cartoon and is called the alpha-carbon. 
Side groups are attached to this alpha-cartoon. The 
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alpha -carbon is connected to the carfeoayl carbon which 
is called C. C is connected to the carhonyl oxygen 
(called O) and to the H of the next residue. The side 
group atoms are given names composed of the symbol for 
the element (C, O, S) , a Greek letter (alpha, beta, 
gasaaa, delta, epsilon, a eta and eta) , and perhaps an 
arable numeral if the side group is forked . 

A thermodynamic quantity which describes the 
behavior of systems in dynamic equilibrium. The symbol 
G represents the Gibhs free energy- 

x.;; "eract jon 
The electrostatic attractive interaction is defined 
in terms of Coulomb's law: 

&E - S a % e*/t>*** r afo 
where & E is the energy change that results from bringing 
two charges., a and b initially separated by infinity,, to 
some distance, r a &. S a and z h are the respective number 
of unit charges? e is one unit of charge. when the 
charges a and b are of opposite signs, the coulomb! c 
interaction is attractive. When the positive charge is 
a divalent metal ion, and the negative charges come from 
an array of points in the structure of a protein, then 
the same energy involved in this attractive interaction 
must be overcome in order to unfold the protein. This 
means that the free energy of metal ion binding is added 
to the overall free energy change for unfolding, A G u , 
thus making this latter parameter more positive and 
shifting the equilibrium to the folded state. 
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A structural segment or segments of polypeptide 
chain folded ia mm a way so as- to give the proper 
geometry and electrostatic configuration for binding a 
divalent metal ion. i$ the physical arrangement of 

protein atoms around a bound metal ion, 

Protain,,,iMfeUi^ 

Protein is defined heroin in terms of the G of 
unfolding. The larger M< the »or® the equili- 
brium folded tr~ unfolded is shifted in favor of the 
folded state. One can got an estimate of from the 

midpoint, of the unfolding transition, Tm f measured by 
scanning oalorimetry . Alternatively one may measure the* 
kinetic stability of a protein from the rate of thermal 
inactivation „ 

This invention is directed to methods for the 
redesign of the tertiary structure of proteins in order 
to increase the stability of their correctly folded and 
biologically and/or eatalytiealiy active state. The 
methods of this invention demonstrate how to increase 
the electrostatic interactions between proteins and 
divalent metal cations in order to increase the binding 
constant for metal ion binding. The Increased electro™ 
static attraction between the protein and metal ion as 
measured by an increase in the free energy of 
dissociation, directly translates into an 

increase in the free eneroy of unfolding, ,&<3 U , for the 
protein of interest, thereby yielding enhanced stability 
for the folded protein structure.. 
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The guide for. designing the alteration of the Betel 
ion binding site through the substitution, deletion f or 
insertion of amino moid residues is Coulomb's law: 

E - S a Z h e 2 /D aff r ab (equation 1} 

where AB, measured in Kcal/aole , is the change in energy 
for a system that consists of two point charges, a and 
b, as they are brought together as a function of their 
distance of separation, r ab , 3 a and S b are the 
respective number of unit charges; « is one unit of 
electronic charge (4,8032 X lG~ i£J esu) f and D is the 
dielectric constant. In the case of a divalent jsetai 
cation binding site, s a * 2*. If Z h is a charge of the 
opposite sign then the interaction is attractive and the 
energy, of the system is decreased, But if the sign of 
Z h is positive then the interaction becomes repulsive 
and the energy of the system is increased. It should fee 
noted that for natural systems like proteins, where the 
charges are not point charges, there is a limit to how 
small -rah can fe® c ome before sterie problems result. The 
lower limit of therefore has been found to be in the 
range of 2 to 3ft for low molecular weight model metal 
ion complexes; . 

The special case where the electrostatic 
interaction is between a metal ion and negatively 
charged or dipolar iigands has been thoroughly studied 
and understood through the use of the electrostatic 
•crystal field theory (CFT) first expounded by H. Beth® 
in 1929, This theory treats the interaction between the 
iaeta.1 ion and the Iigands as a purely electrostatic 
problem in which the Iigands are treated as point 
charges (or point dipoies) . This theory did not 
adequately explain the covalent nature of some metal- 
ligand interactions and has been modified into the 
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adjusted crystal field theory (ACFT) or the. ligand field 
theory (LFT) in order to- allow for coval.e-ncy (Cotton & 

Wilkinson. MSlfflcMJsa^a^Big. S&emitry, Xntersciencs 

John Wiley pub,, H.Y. 3d «d. (197a)), However, in the 
oases where the divalent, metal, ion in question is from 
the Sroup XI aafcals, , W$*+, Ca? + , Sr 24 ", and 
than the CFT- does quite well in predicting geometry,, 
affinity {%) , and otfeer physical chemical parameters of 
these metal ion complexes,. The reason for the CPT 
success in these cases is due to the fact that the d 
electron orbitals responsible for the co-valency of 
transition metal ion complexes are exspty for the Group 
XX divalent metal ions , thereby greatly simplifying the 
correlation of experimental results with theory . 

One major source of uncertainty that is encountered 
in electrostatic interactions in proteins is the 
dielectric constant, Xt yaries nen-unirormly depending 
on shielding and charge effects. Many attempts have 
teen made to simulate the effective screening 
(dielectric value) of electrostatic charges inside the 
tertiary structures of folded proteins (Matthew, Ann. 
Raw ftibphys, Cham. f Ms 387-417 {1085}) * Even the most 
sophisticated models, however, have a 5 keel/mole error 
range (more than 3 pXa units) when attempting to 
calculate intrinsic* pKa values of ion i sable groups 

inside proteins (Russell & WarsheX, J. :Mol Mol,, f 

li§. J389-404 ( 19 S 3 ) ) . Heverthal ess , electrostatic 
interactions are expected to be significant over 
relatively long distances even with moderate screening, 
and vary inversely with the linear distance. 

The electrostatic forces described above are also 
involved In interactions between uncharged but polar 
molecules, but the energy of the interaction is more 
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camplex than that between siapi© ions. The energy 
expression for such interactions also generally varies 
inversely with the distance between such ®ol«cul«s, 
raised to a power usually greater than 1 tout less than 6 
C cr e igh ton , j&s&filBfl-ii ^js&d^i^ .,.m A ^.Moieo.;alaj: 

Properties, Freeman fi Co* , JT.X. {1984}}. This is due to 
the polarissamility of dipoles when in the vicinity of 
electric fields, i.e., Group XX petals, M 2 *, It is thus 
easier to predict results with charged ions. For these 
reasons, it is preferred to use Group II metals, K 2 *, in 
this invention and charged side groups of amino acids 
residues? COO" of Asp and Qlm HH 3 -* of Lys and the 
amimo-tarminus , and the guanidinium cationic group of 
Arg, 

according to the methods of this invention, a -metal 
ion binding site in a protein can foe enhanced to 
increase the affinity for divalent metal ions. In 
general, the methods of this invention comprise the 
foi lowing steps t 

A.) Correlating the metal „ ion bindiaaLjBdaa^ 

M, sjte with some p a ramet er jBtL_jags&&ia stafoXUty." 

This parameter can be the rats of thermal inactivatisn 
derived from kinetic experiments at elevated 
temperatures and varying concentrations of .metal ion. 
Alternatively , one may use thermodynamic measurements of 
the melting temperature, Tsa, derived from calorimetry or 
some other physical assessment of the unfolding transi- 
tion. These measurements not only allow an estimate of 
the degree to which one is successful in affecting 
stability in a desirable way hut also allow an estimate 
of the affinity of this site for divalent metal ions, 

In identifying the saetal ion binding site for 
correlation with a parameter of protein stability one 
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»ay analyse the three dimensional structure of the 
protein of interest or some evolutionaily related 
variant, The thres~d imsns iona I structure may be 
obtained from published sources or determined by fcnown 
methods of X~r»y crystallography , other methods of 
obtaining structural information may Include circular 
diohroism, light scattering * measuring the absorption 
and emission of radiant energy, neutron diffraction, and 

a&&a$L a^&.XSg&Sa&J^ so 

$M^A*imo.~&QiM$, and tfa^^me^l^m^ia^J^teased - These 
amino acid changes may be introduced into the cloned, 
sequenced gene for the protein by the technique of 
oligonucleotide site-directed in vitro mutagenesis* 
This method is described in detail in Bryan 
^^^atji^cad, Scj^JIM (1986), incor- 

porated herein by reference* 

The criteria for selection of amino acids in close 
proximity to the boned metal ion for substitution, 
insertion, or deletion via site-directed mutagenesis Is 
as follows?:: 

Select amino acids for substitution, 
insertion, or deletion as close as possible to the hound 
metal ion, but without introducing steric hinde.ra.nce. 
Amino acids are selected for substitution, insertion, or 
deletion that will optimise the distance and geometry of 
the electrostatic attractive forces at the binding site. 

(a) in substituting amino acids, 
preference should be given to changing those assin© acids 
that are not conserved in evoiutionaily related 
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homologons proteins. Conserved amino acids in 
evolutionally related homologous proteins will generally 
indicate that these conserved asa.no acids are favorable 
to the protein and should be maintained. In contra- 
distinction, amino acid substitutions should be evalua™ 
ted that differ (variable positions) in the area of the 
metal ion binding site hetween the protein to foe altered 
and in evolutional!? related homologous proteins. Thus, 
for example, if the protein to bo altered contains a 
neutral residue in the region of binding and a related 
protein contains an &sp or Glu (negatively charged) , 
then high priority should be given to substituting the 
neutral residue with asp or Glu. 

(b) In insertion and/or deletion of 
amino acid residues, preference should be given to those 
amino acids that will .optimise the distance and geometry 
of ligand interactions . In general,, this involves 
creating a radius between the amino acid(s) end the 
metal, ion as close as possible to 2*$ A that can be 
obtained without creating steric hindrance. For 
example; if a position when changed to asp or Glu is not 
at the optimum distance (2.3 A) from the metal ion, then 
the introduction of an insertion or deletion in this 
vicinity may allow the Asp or Glu to swing in closer in 
distance to that found to be optimal. 

(c) Positively charged residues that are 
close to the bound metal ion should be changed to 
neutral or negatively charged residues. Positively 
charged residues may also simply be deleted. Neutral or 
negatively charged residues may be: inserted. 
Combinations of insertions., deletions, and substitutions 
are contemplated according to the methods of this 
invention. As the divalent metal ion has a positive 
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charge, changing positively charged residues to neutral 
or negatively charged residues will. increase the 
affinity of the metal ion binding site. 

Candidate amino acid substitutions selected fey 
the criteria given above may then fee simulated on a high 
resolution computer graphics system such as the Evans 
Sutherland Model PS330 interfaced with a computer of 
suitable configuration such as the VAX 11/780, such a 
graphics system allows one to graphically make the 
proposed change and test whether there are any obvious 
staric problems that the new amino acid residue will 
introduce f assuming there is no movement of the main 
chain atoms of the polypeptide fold. The criteria 
employed in such an analysis are as follows: 

a) Measure the interatomic distances 
between most of the atoms of the proposed residue's side 
chain with those of most of the atoms of its nearest 
neighbors in the protein structure, 

&•} If there are any serious violations 
of the van der Waais radii of any of the atoms or severe 
electrostatic repulsion between any of the negative 
charges f then rotate one or more of the bonds in an 
attempt to reduce these effects. If either a) or fej 
gives a structure with no serious van der Waals radii 
violations or electrostatic repulsions,- then introduce 
the proposed change into the protein via site--dirscted 
mutagenesis,. If serious van der Waals radii violations 
or electrostatic repulsions persist, then give the 
change a low priority for experimental evaluation* 

Proteins that may foe redesigned according to the 
methods of this invention include those that have M 2 * 
metal binding sites. The majority of proteins, with 
divalent metal, ion binding sites are enzymes with 
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calcium ion binding sites, although- enzymes with other 
divalent metal ibm binding sites are Known. Selection 
of the proteins to be redesigned will, depend, of course, 
on the ultimate contemplated use of the redesigned 
protein, Further, as psd herein the use of the term 
protein is intended to el so mean protein fragments, such 
as polypeptides . 

The largest class of naturally occurring proteins 
Is made up of enzymes* Each enzyme generally catalyzes 
a different kind of chemical reaction., and is usually 
highly specific in its function, Although there may he 
slight variations in a distinct type of naturally 
occurring enzyme within a given species of organism, 
enzymes of a specific type produced by organisms of the 
sasse species generally are substantially identical fit}* 
respect to substrate specificity, thermal stability, 
activity levels under various conditions {e.g* f 
temperature and $>K) , oxidation stability ,, and the like. 
Such characteristics of a naturally occurring or * s Vild~ 
type* enzyme are not necessarily optimized for utiliza- 
tion outside of the natural environment of the enzyme. 
Thus,, it may he desirable to alter a natural 
characteristic . of an enayme to optimise a certain 
property of the enzyme for a specific use, or for use in 
a specific environment according to this invention, 

Preferably, the methods according to this invention 
are used to mutate serine proteases to enhance certain 
characteristics, particularly thermal stability. k 
protease is a catalyst for the cleavage of peptide 
bonds, A serine protease is an enzyme which has an 
essential serine residue at the active site and 
catalyzes the hydrolysis of peptide bonds . serine 
proteases can be inhibited by phony imsthane- 
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sulfonyl fluoride and by diisopropylilu- hate, A 

sxjtbti.li.siB is a serine protaa.se produced by Gram, 
positive bacterid or fungi. These serine proteases may 
include, but are not limited to, those from Bacillus 
species such as subtil lain BPH* from IlsHI^s. 
mMiiaLfeiiDi, subtilisin Carlsberg from BaaUIaa 
.11 che.nl form is, sm&t ills in DY from Bacillus, m, 
sub tills in amylosachari t ieus from Bacil.lus 
am ylose chariticus r and mesenteric s 

proteases » such as protease X, theraomycoiasa, and 
thermitase from ^erraoacfclnoiavce* yyi<R$<£$ may also be 
enhanced in accordance with the invention, as veil as 
mammal ian proteases produced in a bacterial host. 

Other preferred proteins to be redesigned according 
to the methods of this invention include those proteins 
having weak divalent metal ion binding sites ♦ Examples 
of such proteins include alpha amylase, glucose 
isomerase, the^olysin, and neutral protease, generally 
in wild-type proteine, the divalent metal is calcium* 
although binding sites with other divalent metals are 

known, (Craighton, gia&aiaiSi £&GI£&I£<I M>d,,,,Mqle,colar 

£X£J&£&iS£, Freeman & Co., N.Y. (1984)). 

The amino acid substitutions., deletions, and 
insertions may he accomplished by site-directed {i.e., 
oligonucleotide) mutagenesis. Site-directed mutagenesis 
is .well-known in the art and is described in detail in 

Bryan e^_ai . , SssSh^SSlSi.ll,. MM^L^sLt^MA S3; 3743-3743 

(19S6) , incoporated herein by reference. 

Further, in another embodiment of this invention., 
regional directed in vitro random mutagenesis is -used to 
alter the amino acid residue (s) in the loops responsible 
for metal ion binding in the protein of interest, the 
preferred procedure for in_jat.ro. random mutagenesis is 
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described in detail in United States patent application 
serial number 828 f 545 and PCT application 
PdT/XJS S ?/ Q O 3 4 8 , footl* Incorporated herein fey reference. 
The new variants generated by tills procedure can be 
screened or selected for thoee variants that exhibit the 
desired parameter of enhanced protein stability, These 
screened of selected variants can then be further 
characterised first through genetic manipulations to 
identify the change (s) in amino acid residue (a) . In 
addition, these screened or selected variants can be 
further identified by kinetic and thermodynamic 
characterization of stability. Finally, these variant 
proteins may be crystallised and X-ray crystal 
structures obtained at high resolution to fully 
correlate the increased stability with known protein 
structural information. 

To further enhance the stability of these variant 
proteins, additional modifications to the amino acid 
residue (s) in close proximity to a metal ion binding 
site(s) may be effected according to the previously 
described methods of this invention. 

in general, prior to mutation of a gene coding for 
an enssyme of interest, the gene generally is first 
isolated froia its natural source and cloned in a cloning 
vector. Alternatively, mRNA which is transcribed from 
the gene of interest can be isolated from the source 
cell and converted into cDHA by reverse transcription 
for insertion into a cloning vector. A cloning vector 
can foe a phage or piasmid, and generally includes a 
replicon for autonomous replication of the vector in a 
microorganism independent of the genome of the micro- 
organism, k cloning vector advantageously includes one 
or more phenotypic markers, such as m& coding for 
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antifeiotic resistance, to aid in selection of 
microorganisms transformed by the vector. 

Procedures for insertion of SKA or cDHA into a 
vector for cloning purposes are mil kn©*m in the art. 
These procedures generally include insertion of the gene 
of interest into an opened restriction «ndox*ucX«as'e site 
in the vector, end may involve addition of polymeric 
tails of daorynucleotides to the ends of the gene and 
linking the gene to opened ends of a cloning vector 
having complementary p 

& gene of interest to he altered fey site-directed 
mutagenesis can he present in an expression vector, An 
expression vector generally falls under the definition 
of a cloning vector since an expression' vector usually 
includes the components of a typical cloning vector, 
namely one or snore replicons as defined amove, and one 
or asore phsnotypie markers for selection purposes. 
Additionally , an expression vector includes control 
sequences encoding a promoter, operator, ribosome 
binding site, and translation initiation signal. For 
expression under the direction of the control sequences, 
a target gene to fee treated according to the invention 
is operahly linked with the control sequences in the 
proper reading frame. An expression vector containing 
the 0*fA sequence to be targeted can he a phage or a 
plasmid, with piasmids being preferred. 

By using the methods of this invention, subtiiisin 
7172 with an amino acid substitution of Pro 172 to Asp 
has been created. As shown in Figure 4, the new 
subtilisin 7172 has a five-fold increase in the affinity 
for calcium over wild-type sufetilisin BJ?H 5 . Sufotilisin 
S312 was. also created,, with an amino acid substitution 
of Pro i?2 to Gin. Subtil isin S3 12 was found to foe 
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similar to subtilisin 7172 ♦ subtilisin 7148 with a 
substitution of Gly 131 to Asp has also been created* 
Figure 5 shows that hew subtilisin 7X48 has a two-fold 
.increase- in the affinity for calcium over the wild-type 
subtilisin 8PH J . As discussed in the Examples , the dif- 
ference in the binding strength is believed to foe due to 
the distances .between the substitutions and the binding 
sits, in ?r?2 and 8312 it is si and in 714$ it is isA 
away from the divalent metal cation ; 

This invention is intended to be useful for the 
stabilization of many different kinds of proteins under 
a variety of conditions. For example, microbial serine 
proteases are used as an active ingredient in detergent 
formulations for the enhancement of detergent 
performance. Sunt ilis.in~2 ike proteases are currently 
used in detergent formulations to improve the removal of 
protein stains such as blood and milk from soiled 
garments. These detergent f ovulations can often 
present conditions of pH, temperature, free metal ion 
concentrations, and detergent content < hydrophobic ity) 
that do not favor the properly folded state of enzymes,, 
i.e. subtil isin-li3ce proteases. Therefore, the present 
invention provides subtilisin enrymes of enhanced 
Stability that may be used in these applications. 

In particular, subtilisin enzymes of this invention 
can be used as an additive to washing preparations, such 
as liquid detergents, especially citrate-based liquid 
detergents, which are used for cleaning fabric. The 
subtilisin enzymes of this invention are more thermally 
stable than wild-type subtilisin and thus does not lose 
activity as rapidly as wild -type when stored in solution 
with detergents or when subjected to high heat during 
use in cleaning. By use of the subtilisin enzymes of 
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this invention as an additive in washing preparations,, 
the removal of proteinaeeoue stains on fabric is 
improved. The amount of subtil.is.in enzyme that Bay bo 
nned as an additive to washing preparations are wall- 
known in the art, or may readily be ascertained by 
routine e5qp.erimentat.i0n. as optimal range of enzyme 
concentration will., of. warse, be related to the cost 
of the ensyme and the amount of cleaning needed. 
Typically, the amount of subtiiisin ensyme added to a 
washingn preparation will be from about 2000 to about 
4000 Alkaline Delft Units/gram (ABG/gm) of washing 
preparation « 

The following examples are presented to illustrate 
the invention and are not intended to be limiting, 

EXAMPLE I 

OF SUBTXLXSlBf BPM S 

An X-ray cryat&il ©graphic structure of subtlixsin 
BPN* at 1.3A resolution is shown in Figure 1. Two metal 
ion binding sties were found at opposite ends ox the 
globular protein. One site is located near the m$~ 
terminus and is labeled calcium site A since this site 
was found to preferentially bind calcium ions even in 
the presence of EBTA at 4*C h more detailed view of this 
site is shown in Figure 2. Seven prcte in-derived oxygen 
ligands are arranged about site A in an approximately 
octahedral geometry if one considers the two Asp 41 
oxygen ligands as one site. Three of these ligands are 
main chain carhonyl oxygen* atoms that originate from lie 
79 , Val 81, and Leu 75, The other four originate from 
the side chain atoms of Sop 41, Gin 2, and Asn 77* .All 
sever; have an average oxygen -calcium distance of 2,4A, 
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The second metal .binding site is labeled calcium 
site 8 and is .located about 2 OA from site A, The second 
site is located in a crevice defined by two segments of 
the polypeptide chain. The first segment contributes 
three main chain carbonyl oxygen atoms as Uganda; Sly 
169, Tyr 171, m& Val 174 (see Figure 3), and the second 
segment contributes an additional main chain oxygen atom 
(Glu 195) along vitii an oxygen atojs derived from a side 
chain carboxylata group (Asp 19?) . Two ^ater molecules 
round out the seven member coordination sphere. Because 
this site is located on the surface it is more exposed 
to solvent than is site A, The mean cation-oxygen 
ligand distance is also longer in site B (3,9A vs. 2.4A 
in site A) . See Table 11 for eettaai calcium*li$an4 
distances in the two sites. 
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fable II. Cation-Ligand* Distances from 13A 
Resolution X-ray Model , 
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*Digand is defined as those atoms that are in 
the. first coordination sphere of th«t metal 
ion, and which have formal negative or partial 
nefative charges. 

The chemical properties of these ' two sites are as 
different as their structures might suggest* For 
example, a single crystal of wild type sxibt.il is in was 
soaked for 24 hours in 5 jbM EDTA near 4*C and 
diffraction data was collected to 2.QA resolution » 
Difference Fourier [ F (native) -F(EDTA] analysis revealed 
that no change had occurred at site A. This result is 
consistent with metal ion analysis of the wild type 
protein after exhaustive dialysis against a 1 mH 
solution of BDTA (2D sM Hepes) at 4 5 C, Under these 
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conditions it was found that the protein still bound 
0.45 moles of calcium per sole of ensyme, Results 
similar to this had hmm reported toy Voordouw sSlJlL^ 
Biochemistry, aSf371«-37 {1976} who estimated the % of 
calcium binding to subtil isin SPN ! to be less than that 
for EDT& {% 10" 11 } at 25*€, Wo have estimated the % 
for calcium binding to this site to be between 10~ 8 and 
ID" 12 near €Q a C, as judged by differential scanning 
calorimetry in the presence of different metal ion 
chelates. 

The chemical properties of the calcium site B are 
somewhat dif ferent, The F<native>~F(SDTA> - difference 
Fourier analysis indicated a rearrangement of the ion 
environment at site B after treatment of the crystal 
with 5«M B0T3U Interpretation of this difference 
electron density together with observations made with 
other subtilisin crystal preparations have lead to the 
suggestion that the calcium ions at Site B can be 
readily replaced by smaller monovalent cations such as 
sodium. This replacement is accompanied by a small 
shift in ion position and a reorganisation of solvent 
about the ion so that the coordination number is reduced 
from 7 to 5. Binding of other monovalent cations (X*, 
and T1+) has also been reported by Drenth et,,,,al A , 
Moffhem, 26.5l??~l8X (1972)* 

Evidence for a weak calcium binding site first came 
from thermal inaetivation kinetic experiments conducted 
as a function of free calcium ion concentration. As the 
calcium concentration was increased above 0,1 m, the 
half-life for the rate of thermal inactivation at 65'C 
was observed to increase from about 1 minute to 
approximately 230 minutes at calcium ion concentrations 
of 100 m to 300 mM (see Figure 4} . This data was found 
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to closely fit a theoretical ti~rd':ioa curve for a 
binding site that has a Kd of 32 rsH for a single calcium 
ion. Other Group ix metal ions such as Mg 2 *, Ba 24 \ and 
Sr 2 * were similarly tested* The former two metal ions 
wore found to fee ineffective for the stabilisation of 
subtil isin i& this same concentration range as shown in 
Figure 4. The Sr 2 * ions, however f were £©y-«d to 
increase the thermal stabilitf of the emgyma but the 
higher concentration necessary to have an effect 
suggested a much weaker affinity for this site* The 
concentration of SrCl2 necessary to give a stabilisation 
equal to that of CaCl 2 ^ about 20 times higher, 
suggesting a Kd for binding of Sr 2 * to toe in the 0.7 K 
range « These results with the Group II metal ions are 
typical of calcium binding sites in proteins (Stuart 
al«, Mature $24: 84*87 {1986}} and reflects a selectivity 
based on eisse preference. The ionic radius of Mc 2 *, 
0.78 A (hadd radius? Cotton & Wilkinson, MsaaSftSl 
irjoxgani ^Chemisto: , Intersciemoe , John Wiley, N.t* 3d 
ad. (1972)), is apparently too small, and the ionic 
radius of Ba 2 *, is apparently too large, The 

ionic radius of Sr 2 +, X*27i, can be accommodated, but 
the ionic radius of Ca 2 *, l«0ii, apparently fits the 
bast-. A similar result- was found by Bolguin $&J&JU 

Maflfanu MasfeE§u., ili^s-is (ism)) for purified ca- 

ATPase from sarcoplasmic reticulum. 

Xt is also possible to measure the effect of the 
free -calcium ion concentration on the thermal stability 
of subtil isin BPH* through the use of differential 
scanning .calorimetry (DSC) . Increasing the calcium ion 
concentration above 0,1 mM causes the midpoint for the 
thermally induced unfolding transition, Tm, to increase 
from about m,S v C for the protein as isolated where it 
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binds 1 mole calcium/tsola protein, to approximately 8(TG 
for the protein in the presence of 150 mK Ca€l 2 . 

The next necessary step taken was to correlate the 
data in Figure 4 with one of the two structural sites 
identified by X-ray cry stenography . In light of all 
the available chemical and physical information 
described above, it was possible to nominate the calcium 
site 8 as the most likely candidate responsible for the 
weak calcium binding demonstrated in Figure & for wild 
type sufctilioin BFN* « 

Through the use of the rules outlined above for the 
selection of the amino acid residues to be altered, a 
decision was made to change Fro 172 to Asp, Pro 172 is 
part of one of the loops that comprise the calcium site 
3 , The simulated earboxylate chain is approximately 
5A from the fecund metal ion, this particular amino acid 
residue is somewhat variable in related subtil is in 
sequences. (It exists as an Asp in subtil is In 
Carisfoerg-) finally, simulation on a graphics screen 
showed that it should fit without van der Waal 8 radii 
violations. The change was introduced by site-directed 
mutagenesis t essentially as described above and by Bryan 
et a"l., P**c. Nat»l Acad. Sci. USA, 3743-3745 (1986) 
(incorporated herein by reference), and the new strain 
of S,. c -:^:y.,ls containing the altered subtil lain gene 
was called 0X717 2. The variant protein product of that 
gene was purified and analyzed for enhanced metal ion 
binding and stability. The results are shown in Figure 
4 and compared with that for the wild type protein (the 
new protein is labeled 7172), 

The calcium ion dependence on the half-life (tl/2) 
for thermal inaetivation of suhtilisin 7172 is clearly 
shifted toward lower calcium Ion concentrations. The 
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theoretical titration csorvre that best approximates the 
data is one for a single sit© with % » € m* Xt 
therefore appears that changing Fro 172 to Asp increases 
the affinity of the B site for calcium by about 5 ~ fold 
over the wild type.. The mgnitude of this change, ^pK «* 
0*7 , is consistent with theoretical considerations for 
the introduction of a negative charge at about si from a 
metal cation given the known difficulties currently 
associated with the application of electrostatic theory 
in proteins. Further support for this supposition comes 
from mutagenesis experiments with staphylococcal 
nuclease where it was found that changing the Asp 40 
earboxyiate ligand to &ly caused an 8 -fold decrease: in 
calcium binding to the weak calcium sits (Kd * 1 «H) in 
this protein (Serpersu et^l. , Biochemistry 26:1283-1300 
<19S?})*- Those results therefore provide convincing 
evidence that the weak calcium biding site measured toy 
thermal inactivation corresponds to the calciuja 8 site 
identified through X-ray crystallography* Moreover, 
these results clearly demonstrate the anility to 
Stabilise subtil is in BPH ? in the presence of decreasing 
concentrations of calcium by increasing the the metal 
binding affinity of the site. The extrapolated outcome 
of continued alteration of site S is the eventual shift 
in the titration curve so that % is so small that no 
excess metal ion {or dramatically reduced amounts}- will 
be necessary to produce the stability one sees at 100 a*M 
Ca for the wild type protein. 

In addition to changing Fro 172 to Asp, Fro 171 was 
also changed to Gin- The calcium dependence of the 
rates of thermal inactivation were also examined for 
this new variant {8312} and the results were found to be 
almost identical with that of subtilisim 7172, 
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The amino acids at the calcium A were also 

changed to increase thermal stability. In strain 8347, 
As« 7S was changed to Asp. In strain 8364, Ser 78 was 
changed to Asp. In strain 8374, these mutations, Asu 76 
to Asp and Ser 78 to Asp, were, combined by 
oiigonucleotide-directed mutagenesis. As shown in Table 
1X1, the calcium binding affinity of 8374 is increased 
in an additive fashion over each mutation alone. All 
three strains showed increased stability over the wild- 
type (BPS*) and over suhtilisin Carlsberg. 

Furthermore, the changes made in the calcium B 
binding site may be made in combination with those made 
in the calciins A binding site in order to increase 
stability. Thus, for example, the mutation Asn 76 to 
Asp and Pro 172 to Asp (or Glu) may be made. Additional 
examples would include Asn 76 to Asp, Ser 78 to Asp, Pro 
172 to Asp (or Glu) and Gly 131 to Asp? Asn 76 to Asp 
and Gly 131 to Asp? Ser 78 to Asp and Gly 131 to Asp ; 
Ser 78 to Asp and Pro 172 to Asp (or Gin } t and ser 78 to 
Asp, Gly 131 to Asp and Pro 172 to Asp for Glu) , 
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EXAMPLE II 

Insertions and/or Deletions of Amino Acid Kesidtt.es- • 

In addition to substitution of ami.no acid 
residues, insertions and/or deletions of ani.no acid 
residues in the vicinity of the binding site of interest 
can be done in a tether attempt to influence the 
electrostatic forces at that site. This can be thought 
of as a method for optimising the distance and geometry 
of iigand interactions. For example,, an insertion 
between residue 172 and 173 in one of the loops that 
comprises the calcium B site in subtiiisin. 7X72 may he 
effected in order to move the Asp at residue 172 closer 
to the bound metal so that it may become a formal Iigand 
and possibly displace one of the bound water molecules. 

In searching the Brookhaven Protein Data • Bank for 
peptide sequences the con forma ti onal structure of which 
match parts of the subtil is in BPJsT loop from residue 168 
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to residue 175, several sequences were found which 
matched well at residues 170-171 and at residues 173- 
174, but ■which haa. one extra aiaino acid residue 
inserted. Ail of these peptides had the main chain 
deformed in a way that would permit the side chain of 
residue 172 to tip mors toward the calcium binding 
sit*. The carbonyl oxygen* which bind to calcium were 
not significantly disturbed in this computer simulation. 
If the residue inserted happens to be an asp or a Gin 
then one would gain another opportunity for further 
influencing the electrostatic interaction at this site. 

An amino acid deletion may also cause a mora 
favorable environment for metal ion binding by altering 
the intercharge distance and geometry of neighboring 
residues, 

EXAMPLE 111 

Enhancing Stability Using Random Mutagenesis 
The 7142 variant of subtil is in was generated and 
identified by a .method of random mutagens is which 
exposed the entire subtil isin gene to mutagenic agents » 
(See United States Patent Application serial number 
828,545 and POT Application PCT/US87/0034-8 . J In this 
particular case the mutagenic agent was sodium bisulfite 
and the filter assay identified this protein as a stable 
mutant. DNA sequencing revealed that subtil isin 7142 
had two amino acid changes: Sly 131 to Asp, and Ala 116 
to Thr. These changes were then introduced individually 
by site-directed mutagenesis in order to investigate the 
physical properties of the individual amino acid 
substitutions, the Sly 131 to Asp change alone was 
found to be solely responsible for the observed stabil- 
ization. This new variant was designated subtil is in 
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7148, ThB Ala 11$ to Thr change was found to have no 
effect on the rate of thermal inactivation. By changing 
Sly 131 -to Glu, the results on calciuas ! . ien >s of the 
rates of thermal motivation may be similar with the 
amino acid substitution to tec* 

When suhtilisin 7X4B was analysed for Kinetic 
stability as a fuucfcion of the calcium ion eo ncent r a t i on 
it was found to have a slightly increased affinity far 
calcium ions as shown in Figure s« .In this case, the 
data was found to tee best approximated by a theoretical 
titration curve for a single binding site with a Ed of 
15 mM. This apparent a~folu increase in metal ion 
affinity, or &pK « 0,33, correlates very well with X-ray 
crystallography data on subtilisin 7142. The 1*8 t 
resolution crystal structure for subtilisin 7142 shows 
the earbossyiat© group of the new Asp 131 residue to be 
approximately 10 1 from the metal ion at the calcium B 
site. Therefore the magnitude of the observed increase 
in metal ion affinity is about what one would expect in 
light of this intereharga distance. For example, Thomas 
et.al .. Nature US; 375-376 (1985) , observed a 0.29*0.04 
shift in the pKa of His 64 of subtilisin BPH» after 
changing Asp 99 to ser which is 14 1 away. Moreover,, if 
the introduction of a negatively charged residue 5 A 
away from the metal ion bound at the calcium 8 sits 
causae a ApK * 0.7* then it is entirely consistent with 
the theory (equation 1} that the introduction of a 
negatively charged group at 10 i from the site gives a 
iipK «* 0.33. 

The mutations €iy 131 to Asp and Pro 172 to Asp 
were combined by oligonucleotide-directsd mutagenesis to 
create strain 8331. The calcium binding affinity of 
S331 is increased in an additive fashion. 
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Regional-directed In ,v.ifcxg random mutagenesis can 
be used to alter sutotilisin BPB< for enhanced thermal 
stability. In this particular application o.f random 
mutagenesis only that part of ifc* m&till&Ln gene that 
cooes for asi.no acid residues in the loops responsible 
for Metal ion binding is subjected to mutagenesis. The 
methods used for random mutagenesis are those described 
in United States patent application 828,545 and FCT 
application PCT/US87/0Q348 , both incorporated herein by 
reference* 

in the subtil isin B1?H* t the calcium B site is 
comprised of only two short stretches of polypeptide 
chain (residues 169-174 and 193-19?}. Because of this,, 
it is relatively convenient to specifically target this 
region for in vitro random mutagenesis. This was 
accomplished toy creating recognition sites for restric- 
tion endonucleases in the cloned subtil lain gene at 
sequence flanging the two target areas. For example., to 
target the residues 169 to 174 loop of the B-site, 
oligcnucleotide-directed mutagenesis was used to create 
an XhoX cutting site at residue 162 and a SalX site at 
residue 176, By taking advantage of the redundancy of 
the genetic code, these restriction sites were 
introduced without altering the amino acid sequence 
specified by the gene. Random mutations were created in 
the target region toy cleaving the gene with lhp.1 and 
Sgll and then ligating synthetic DKA fragments into the 
gap created in the gene by the two restriction ensymes. 
The DMA fragments were synthesized in such a manner that 
changes from the native DMA sequence occurred at random. 
Using this procedure, the synthesis of the DHA insert 
can be controlled so that mutations occur at any 
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frequ&ncy, from wry minor variation from the native 
sequence to a completely random sequence. 

Although the foregoing invention has hesm described 
in some detail by way of illustration and example for 
purposes of clarity and understanding, it will foe 
obvious that certain changes and modifications may he 
practiced, within the scope of the invention,, as limited 
only by the scope of the appended claims* 
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1, h method for increasing the stability of a 
protein containing a metal ion hin$im comprising 
the following steps? 

{1} correlating said metal ion binding site 
of said protein with a parameter of protein stability 
ana 

(2) altering the structure of said metal ion 
binding sit* through the substitution, insertion, and/or 
deletion of assi.no acid residues in close proximity to 
said site so that the electrostatic attraction between, 
amino acid residues at the mt&X ion binding site and 
the metal ion is increased. 

2, The method of claim 1 wherein said correlation 
of the metal ion binding site with a parameter of 
protein stability is accomplished by identification of a 
metal ion binding site through analysis of the three- 
dimensional structure of said protein or an 
evolutional ly related homologous variant. 

3, The method of claim 1 wherein said altering is 

accomplished through oligonucleotide-directed in yjtj^ 

mut agenesis « 

.4. The method of claim 1 wherein said altering is 
accomplished by in vitro random mutagenesis and then 
screening or selection of variants that exhibit a 
desired parameter of enhanced protein stability. 
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The method of cleim 1 wherein said step (3) 
further comprises: 

examining asdno acids of said metal ion 
binding structural sit* aM selecting amino acids for 
substitution., insertion, or deletion that are as close 
as possible to the hound metal ion, without creating 
unacceptable stetio hindrance, and that optimises the 
distance and geometry: of the electrostatic attractive 
forces at the metal ion binding site, 

6* The method of claim 3 wherein the amino acids 
selected for substitution, insertion,, or deletion are 
ones which are not conserved in evolutional ly related 
horn o i ogous pr ote i ns » 

7, The method of claim 5 wherein neutral or 
negatively charged amino acids are substituted for 
positively charged amino acids. 

8, The method of claim 3 wherein neutral or 
negatively charged amino acids are inserted into the 
metal ion binding site. 

9, The method of claim 5 therein positively 
charged amino acids are. deleted from the metal ion 
binding site. 

10, & natal ion binding site-containing protein 
with increased stability over the wild-type protein 
obtained by altering the amine acids at the metal ion 
binding site of said protein to increase the electro- 
static attractive interactions between the amino acids 
at said metal ion binding site and the metal ion. 
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IX, a jsetal ion binding site-containing protein 
with increased stability ovar the wild-type protein 
produced according to the method of claims 1 to §, 

12. The protein of claim 10 comprising an ensyma, 

13. The protein of claisi 12 comprising a serine, 
protease. 

14. The protein of claim ID wherein said metal ion 
binding site is a divalent metal ion binding site and 
wherein the metal ion is selected frois ca 24 \. Hg 2 +, Ba 2 *, 

15. The protein of claim 14 wherein said metal ion 
is Ca 24 l 

16. A serine protease with increased stability 
over the wild-type serine protease obtained by altering 
one or sore amino acids at the calcium ion binding site 
to increase the electrostatic attractive interactions 
between the atsin© acids at the calcium ion binding site 
and the calcium ion. 

17. The serine protease of claim 16 comprising 
snbtiliein, 

18 x The serine protease of claim 17 wherein said 
snbtilisln is selected from the group of wild-type 
subtilisin consisting of subtilisin BPS', subtilisin 
carlsberg, subtilisin W, subtilisin amylosaoharitricus, 
and mesentericopeptidase. 
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19. h subtil is in with increases stability over the 
wild~typ® suhtiiisin obtain by altering the amino acids 
at the calcium binding sits to increase the 
electrostatic attractive interactions between the amino 
acids at th« calcium ion binding site and the calcium 
ion. 

20. h subtil Isin vita increased stability over the 
wild-type protein produced according to the method ©f 
claims % to S. 

21. The- sufotilisin of slain 14 wherein said 
subtil isin has one or more amino acid substitutions 
selected from the group consisting of Aen 76 to Asp , Ser 
78 to mp, Sly 131 to Asp f and Pro 17a to Asp or <*lu. 

22. The serins protease of claim is wherein said 
subtiiisin is coded for by a mutated subtil isin gene, 

23. Subtilisin 7172, with an amino acid substitu- 
tion of Pro 172 to Asp. 

24. Subtilisin 7148, with an amino acid substitu- 
tion of Gly 131 to Asp. 

25. Substilisin 3312, with an amino acid substitu- 
tion of Pro 172 to Glu« 

26. Subtilisin S331, with amino acid substitution 
of Sly 131 to Asp and Pro 172 to Asp. 
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27. subtil iaia S34? f with amino acid substitution 
of jam 76 to Asp, 

28. Subtilisin with amino acid substitution 
of ser 7$ to Asp, 

29.. Subtilisin 8374, with amino acid substitution 
of Asn 76 to Asp and Ser 7S to Asp. 

30. A w*t*% ion binding site-containing protein 
with increased stability over the wild-type metal ion 
binding site containing protein wherein the amino acid 
sequence of the metal ion binding site of said metal ion 
binding site-containing protein differs by at least one 
amino acid from the as in© acid sequence of the metal ion 
binding site of said wild-type metal ion binding site 
containing protein > 

31. The metal ion binding site-containing protein 
of claim 30 wherein the difference in the amino acid 
sequence comprises a substitution of an amino acid. 

32. The metal ion binding site-containing protein 
of claim 30 wherein the difference in the amino - acid 
Sequence comprises a deletion of m amino acid. 

33. The metal ion binding site-containing protein 
of claim 30 wherein the difference in the amino acid 
sequence comprises an addition of an amino acid. 

34. A washing preparation comprising the 
subtilisin of any of claims 17-29 at a concentration of 
2>000 to 4,000 Alkaline Delft Units. 
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35, A method for improving the removal of 
protainaceous stains on fabric comprising contacting 
said fabric with th« washing preparation of claim 34 and 
cleaning said st&lnscl fabric with said washing- 
preparation. 
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